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ABSTRACT 

A analysis of standard model accretion disk synthetic spectrum fits to com- 
bined FU SE and STIS spectra of V3885 Sagittarius, on an absolute flux basis, 
selects a model that accurately represents the observed SED. Calculation of the 
synthetic spectrum requires the following system parameters. The cataclysmic 
variable secondary star period-mass relation calibrated by Knigge in 2007 sets 
the secondary component mass. A mean white dwarf (WD) mass from the same 
study, that is consistent with an observationally-determined mass ratio, sets the 
adopted WD mass of O.TMq, and the WD radius follows from standard theoret- 
ical models. The adopted inclination, i = 65°, is a literature consensus, and is 



^ Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space 
Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, 
Inc. under NASA contract NAS5-26555, and the NASA-CNES-CSA Far Ultraviolet Explorer, which is 
operated for NASA by the Johns Hopkins University under NASA contract NAS5-32985 
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subsequently supported by analysis. The mass transfer rate is the remaining 
parameter to set the accretion disk Tgfj profile, and the Hipparcos parallax con- 
strains that parameter to M = 5.0±2.0xl0~^ Mq yr~^ by a comparison with 
observed spectra. The fit to the observed spectra adopts the contribution of a 
57, 000±5000K WD. The model thus provides realistic constraints on M and Tes 
for a large M system above the period gap. 

Subject headings: Stars: Novae, Cataclysmic Variables, Stars:White Dwarfs, Stars: 
Individual: Constellation Name: V3885 Sagittarius 



1. Introduction 



Cataclysmic variables (CVs) are semi-detached binary stars in which a late rnain se - 
quence star loses mass onto a white dwarf (WD) by Roche lobe overflow (IWarnerl Il995l ). 
In non-magnetic systems the mass transfer stream produces an accretion disk with mass 
transport inward and angular momentum transport outward, driven by viscous processes. 
The accretion disk may extend inward to the WD; the outer boundary extends to a tidal 
cutoff limit imposed by the secondary star in the steady state case. If the mass transfer 
rate is below a certain limit, the accretion disk is unstable and undergoes brightness cycles 
(outbursts), and if above the limit, the accretion disk is stable against outbursts. The latter 
objects (which have no rec orded outburst of any t y pe) are called nova-like (NL) systems. As 



shown bv current st udies (IRibiero fc DiazI (120071): iHartlev et al.l (12005. 120021): IWoods et al. 



(1992): lMetzl (Il989h : iHaug fc Dreschell (Il985f ): ICuinan fc SionI (Il982[ ): ICowlev fc Crampton 



( 1l977l )) V3885 Sgr is a NL system (see W arner 1995 for a detailed discussion). V3885 Sgr 
has a Hivvarco s (iPerrvman et al.l Il997l) parallax of 9.11 ± 1.95 n ias. The orbital period 
is 0.20716071^ jRibiero fc DiazI l2007h a nd the inchnation z~65° dflaug fc Dreschell Il985 
Hartlev et allbood : iRibiero fc Diajl2007h . 



NL systems are of special interest because they are e xpected to have an accretion disk ra- 
dial temperature profile given by an analytic expression (IFrank. King fc Rainelll992l . eq.5.41) 
(hereafter FKR) which defines the so-called standard model; the expression includes the 
mas s transfer rate M as an explicit var iable. In NL systems, which are above the period 
gap (IHowell. Nelson, fc Rappaportll200ll ). the accretion disk dominates the system spectrum. 
Fitting a synthetic spectrum, based on the analytic model, a proxy of the accretion disk tem- 
perature profile, to an observed spectrum potentially determi nes M, a physical parameter of 



basic importance since it controls the evolution of CV systems ( IHowell. Nelson, fc Rappaport 



20011 ). But the analytic expression also is an explicit function of the WD mass, M^d, and it 
must be determined independently. 
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By the basic paradigm of CV evolution (IHowell. Nelson, fc RappaportI 120011) there is a 
close relation between orbital period and donor mass. In a recent paper, iKniggd (120071 ) has 
calibrated the relations among P, M2, R2, T'eff,2; and donor spectral type. Thus, given an 
orbital perod below the tabular upper limit of 6*^, M2 can be determined independently of 
other parameters. If there is an observationally-determined mass ratio, the WD mass follows 
directly. 

A normalizing factor is necessary in comparing synthetic spectra with observed spectra 
and the factor depends on the distance to the system. If the distance is known independently 
an absolute determination of M potentially is possible. 

A complication in fitting a synthetic spectrum to an observed spectrum is the fact that 
the WD may contribute detectably to the system spectrum; evaluating the WD contribution 
requires th e WD T^ff. A lirnited number of CV systems have known values for their WD Tgg 
(jSiorull999l : IWinter fc Sionll2003l ) with extremely few being NLs. Because of the large M in 
NL systems, the accretion disk dominates the system spectrum in the UV and optical, and 
makes determination of the WD T^s difficult. Thus, the T^g of an additional NL WD is of 
obvious interest. 



2. The FUSE and STIS Spectra 



Table 1 present s the observations lo g. The STIS spectra are the same as three STIS 
spectra described by [Hartley et al.l (120021 ) . 



2.1. The FUSE Spectrum 

FUSE is a low-earth orbit satellite, launched in June 1999. Its optical system con- 
sists of four optical telescopes (mirrors), each separately connected to a different Rowland 
spectrograph. The four diffraction gratings of the four Rowland spectrographs produce four 
independent spectra on two microchannel plates (detectors). Two mirrors and two gratings 
are coated with SiC to provide wavelength coverage below 1020A, while the other two mir- 
rors and gratings are coated with Al and LiF. The Al+LiF coating provides about twice the 
reflectivity of SiC at wavelengths >1050A, and very little reflectivity below 1020A . These 
are known as the SiCl, SiC2, LiFl and LiF2 channels producing a total of 8 (partially over- 
lapping) spectral segments: SiCla, SiClb, SiC2a, SiC2b, LiFla, LiFlb, LiF2a, LiF2b for 
each exposure. 

Four FUSE observations of V3885 Sgr prospectively are available for analysis. An 
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observation on Aug. 10, 2002 has a short exposure time and has inferior S/N ratio compared 
to other exposures. An exposure on 24 May 2000 (UT:01:51:55) has the longest exposure 
time and occurs within 25 days of our chosen STIS spectrum. We chose this spectrum for 
analysis. The exposure used the 30"x30" LWRS Large Square Aperture in TIME TAG 
mode. The telescope collected 11 exposures obtained during 11 consecutive orbits of the 
FUSE spacecraft, totahng 12865s of raw exposure time for the 88 spectral segments; we 
downloaded the spectra from the MAST archive and processed them as described below. 



The data were processed with the latest and final version of CalFUSE (v3.2.1) (jPixon et al. 



20071 ). The main change from previous versions of CalFUSE is that the data are maintained 
as a photon list (the intermediate data file - IDF) throughout the pipeline. Bad photons 
are flagged but not discarded, so the user can examine, filter, and combine data without re- 
running the pipeline. Event bursts (short periods during an exposure when high count rates 
are registered on one or more detectors) are processed automatically in this latest version of 
CalFUSE. The bursts exhibit a co mplex pattern on th e detector, their cause, however, still 



is unknown (it has been confirmed (jPixon et al.l 120071 ) that they are not detector effects). 



The 11 exposures all have essentially identical continuum flux levels and there is no 
detectable change in flux from exposure to exposure. The main features that vary from 
exposure to exposure are the depth, width and shape of the absorption lines. The next 
step is to coadd the 11 exposures, for the 8 individual segments. The count rate indicates 
screening of 22s of high voltage (HV) for the 8 spectral segments. In addition there was a 
loss of good exposure time due to event bursts of 472s, 409s, 390s, and 511s in segments 
la's, Ib's, 2a's and 2b's respectively. We obtain a total of 12310s, 12373s, 12392s and 12271s 
for segments la's, Ib's, 2a's and 2b's, respectively, including more than 8200s of NIGHT 
time (i.e. when the spacecraft is located in the shadow of the earth). Since the total good 
collection time varies, we weighted the segments accordingly (see below). 

During the observations. Fine Error Sensor A (which images the LiFl aperture) was 
used to guide the telescope. The spectral regions covered by the spectral channels overlap, 
and these overlap regions are then used to renormalize the spectra in the SiCl (a & b), 
LiF2 (a & b), and SiC2 (a & b) channels to the flux in the LiFl (a & b) channel. The 
low sensitivity portions (usually the edges) of each channel are discarded. In addition, there 
exists a narrow dark stripe of decreased flux in the spectra running in the dispersion direction, 
affectionally identified as the "worm", which can attenuate as much as 50% of the incident 
light in the affected portions of the spectrum, due to shadows of the wires in a grid above 
the detect or. Because of the temporal changes in the strength and position of the worm. 



CalFUSE (IDixon et al.l 120071 ) cannot correct target fluxes for its presence. Therefore, we 



carried out a visual inspection of the FUSE channels to locate the worm and we manually 
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discarded the portion of the spectrum affected by the worm. :~1130-1189A in the LiFlb 
channel. 

We combined the individual exposures and channels to create a time-averaged spectrum, 
weighting the flux in each output datum by the exposure time (as mentioned earlier) and 
sensitivity of the input exposure and channel of origin. The final product is a spectrum that 
covers almost the full FUSE wavelength range 905A-1187A; because we disregarded that 
segment of the spectrum affected by the worm, the spectrum ends at ~ 118 2 A. The absorp- 



tion li nes in the final spectrum are blue shifted by about 0.4A. Based on the lRibiero fc Diaz 



( 120071 ) ephemeris, the FUSE spectrum extends from orbital phase 0.29 to orbital phase 1.01. 



2.2. The STIS spectra 

The STIS spectra (from the MAST archive) were taken in TIME- TAG mode, with the 
FUV-MAMA detector. The E140M echelle grating was used (with a central wavelength of 
1425A), and the 0.2x0.2 aperture. The STIS spectra were processed with CALSTIS version 
2.22. The calibrations were performed using the STSDAS package. Each STIS spectrum 
includes 44 echelle spectra, each one covering roughly 18A. The echelle spectra overlap except 
in the longer wavelengths. As a consequence when the echelle spectra are merged to produce 
the final spectrum, there are 4 gaps in the longer wavelengths. The gaps are about ~0.5A 
wide and are at intervals of about 18A. We have removed the gaps by hand. 

The first and third STIS spectrum have closely similar continua and differ slightly in 
the shape of their absorption lines, with a blue shift of l.SA and 0.6A. The second STIS 
spectrum has a continuum flux level 14 % larger and its absor ption lines are blue shifted 



by only 0.4A. Based on the ephemeris of iRibiero fc Diaa (120071 ). the Table 1 STIS spectra 



extend in orbital phase from 0.67 to 0.81, 0.33 to 0.47, and 0.35 to 0.49 respectively. The 
STIS spectra cover the interval from 1150A to 1720A. 



2.3. The combined spectra 

The FUSE spectrum and the first STIS spectrum of V3885 Sgr are shown in Figure 1 
before correction for reddening; the spectra are characterized by broad and deep absorption 
lines from the source (disk plus wind/chromosphere). These broad absorption lines are blue 
shifted by about 1.5A and 2.5A in the FUSE and STIS spectra respectively; they are about 
4A in width (~l,000km/s) and exhibit the same asymmetry in both the FUSE and STIS 
spectra. The red edge of the absorption is steep, while the blue edge has a more gradual 
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slope, resembling the blue-shifted absorption part of a P-Cygni profile. However, only the 
Civ (1550) and Nv (1240) lines exhibit a true P-Cygni profile, while all the other broad 
absorption lines do not. All these lines belong to high ionization species: Nlll, Niv, Nv, 
Sim, Si IV, Cm, Civ, Ovi, Pv, and Heii. 

In addition to the broad lines, there are many narrow and deep absorption lines. They 
have a width varying between about 0.15A (e.g. Ol) and O.SA (for Hi) and are also blue 
shifted, but by a much smaller amount (~ 0.2A ) than the broad lines. Many of the narrow 
lines belong to low order ionization species; however, some also belong to higher ionization 
species. Some of the sharp lines are known to be interstellar, such as An and Alll. For 
simplicity we put all these lines in t he ISM category. Additional details concerning the STIS 



spectra are in [Hartley et al.l (120021 ). Table 2 and Table 3 list all the lines, separated into 
FU SE lines and STIS lines. 

From the AAVSO web sit^ we determine that during the month of April 2000 and 
through May 2000, the V3885 Sgr visible magnitude varied by only 0.2. The first of the 
Table 1 STIS spectra falls within this time frame and it is closest in time to the FUSE 
spectrum. Also, the orbital phases of the two spectra are nearly the same. Finally, the first 
STIS spectrum most closely matches the FU SE spectrum in the overlap region. All of the 
spectra have been corrected for interstellar reddening of E(B-V)=0.02 (§3). We choose the 
first STIS spectrum to combine with the FUSE spectrum for analysis. Figure 2 shows the 
combination of the FUSE spectrum with the three STIS spectra. The upper panel shows 
the overlap region. The FUSE spectrum is red, the first STIS spectrum is cyan, the second 
STIS spectrum is magenta, and the third STIS spectrum is brown. For optimum fit to the 
FUSE spectrum, the STIS spectra have been divided by 0.98, 1.12, and 0.94 respectively. 
The top panel shows the very good match of the spectra in the overlap region, while the 
lower panel shows that, when normalized, the three STIS spectra track each other with 
excellent fidelity. The STIS spectra are noisy in the overlap region while the upper panel 
of Figure 2 shows that the STIS and FUSE spectra continua superpose accurately. We 
produced a combined spectrum for analysis by terminating the STIS spectrum at 1181. 5A at 
the short wavelength end, dividing it by 0.98, and adding the FUSE spectrum, beginning 
at II8I.5A and extending to shorter wavelengths. 

We suggest, without detailed analysis, that the differences among the STIS spectra 
could be explained by slight variation in the mass transfer rate from the secondary star. 
Whatever the cause of the variation, the fact that the normalized STIS spectra have the 
same SED simplifies calculation of a system model. 



^ www. aavso.org 
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Initially adopted system parameters 



Table 4 (and notes) lists initially adopted parameters and their sources. Several of 
the system parameters are poorly known. As mentioned earlier, t he Hipp a rcos p arallax 
is 9.11±1.95 mas, which corresponds to a distance of 110±30 pc. iKniggd (120071 ) deter- 
mines M2 = O.475M0 for P = 0.20716^^, and a mean M„d = O.7M0. We adopt both 
masses and n ote that the resultin g q = 0.68 is in agreement with t he ob servationally deter- 
mined g~0.7 (IHartley et al.ll2002l ). iPanei. Althaus. fc Benvenutd (120001 ) lists a WD radius 
of 1.091 xlO~^-R0 for a O.7OM0 homogeneous zero temperature Hamada-Salpeter carbon 
model WD. In our subsequent study of the observed spectra we correct the radius for the 
adopted WD T^fi; the contribution of the WD to the system synthetic spectrum is apprecia- 
ble. Tablenotes to Table 4 list alternative sources for the adopted initial paranieters. Those 
notes do not include the determination of i~25° — 30° (ICowley &: CramptonI 119771 ). This 
determination was based on radial velocity observations with a large scatter; the derived 
velocity amplitud e was K = 74 km while much more accurate data gives a value of 
K = 165 km s-i Jflartlev et al.lboosh . 



VerbuntI (119871 ) determined a rough upper limit of E{B — V) of 0.04 (and a preferred 
value of E(B - V) = 0) for V3885 Sgr based on a study of the 2200A "bump" in lUE 
spectra. ISruch fc Engel Jl994[ l hst E{B - V) = 0.02, which we adopt. 



A nurnber o f studies have considered the t idal cutoff boundary, r^, of accretion disks 

( Paczvnski ( 1977 ): Papaloizou fc Pringle ( 1977h : Whitehurst ( 1988 ): Schwarzenberg-Czerny &: Rozyczka 
(Il988h : Iwhitehurst fc Kinej Jl99lh : ICoodmanl (|l993|)). These authors agree on ra^O.SSD, 
where D is the separation of the stellar components. We adopt this expression for the tidal 
cutoff radius of the accretion disk. 



4. The analysis program: BINSYN 



Our analysis uses the program suite BINSYN (ILinnell &: Hubenyl 119961 ): recent papers 
( iLinnell et al.l 120071 . l2008al ) describe its application to CV systems in detail. Briefly, an 
initial calculation produces a set of annulus models for a giv en WD mass, radius, and mass 
trans fer rate. This calculation uses the program TLUSTY|j ( Hubenyl 1 19881 : iHubeny fc L 



anz 



I995I ). The program considers the radial and vertical structure of the disk independently; the 



radial structure is based on the standard model (FKR) and so follows the prescribed relation 
between local T^s and the annulus radius. The vertical structure is solved, self-consistently. 



^ lit t p : / / nova .as tro.umd.edu 
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as described by iHubenyl (Il990l ) and iHubeny &: Hubenyl (Il998l ). The set of annulus models 
covers the accretion disk from its innermost (WD) radius to just short of its outermost (tidal 
cutoff) radius. 



The primary source of vis cosity in CV accretion disks is MRI (IBalbus fc Hawleylll99ll ). 
Hirose. Krolik. fc Stond (120061 ) calculate MHD models with local dissipation of turbulence 
and sh ow that the vertical extent of an annulus is greater than in previous models. iBlaes et al. 
( I2OO6I ) show that magnetic support has a significant effect on synthetic spectra of black 
hole (BH) accretion disk annuli and illustrate the effect in the case of a BH system with 
Mbk = 6.62M0 and with an adopted a = 0.02. They are able to simulate the effect of 
magnetic support within the TLUSTY framework by adjusting the TLUSTY parameters (0 
and (i and use this simulation to calc ulate the effect on syn t hetic spectra. The primary 
effect is a hardening of the spectrum. iKing. Pringle.fc Stond (120071 ) point out that MHD 
models that produce viscosity via MRI require a values that are a factor 10 smaller than the 
a=0. 1-0.4 required by observational evidence and suggest that some caution still is needed 
in accepting the MHD results. Our models have used the default values Co = Ci = 0.0 and 
we assume the effect of magnetic support on our annulus spectra, with M^d = 0.7Mq, is 
insignificant. 

Table 5 lists properties of annuli calculated with TLUSTY(v.202) for a mass transfer 
rate of 5.OxlO~^M0/yr~^. This v alue initially was ch osen for test purposes because of the 
similarity of V3885 Sgr to IX Vel (ILinnell et al.l 120071 ): see §5 for a detailed discussion. The 
control file to calculate a given annulus requires a radius of the WD in units of Rq. All of 
the annuli used the radius of a zero temperature WD for a homog eneous carbon Hamada- 
Salpeter O.7OM0 model from iPanei. Althaus. fc Benvenutd (120001 ). All of the annuli are 
solar composition H-He models, and the models through r/r„d o = 20.00 are converged LTE 
models. The remai ning annuli are LTE - grey models (see the TLUSTY Users Guide for 
an explanation). A IShakura fc SunyaevI (119731 ) viscosity parameter a = 0.1 was used in 
calculating all annuli. Each line in the table represents a separate annulus. Temperatures 
are in Kelvins. The column headed by mo is the column mass, in gm cm~^, above the central 
plane. The columns headed by Zq and iVg are, respectively, the height (cm) above the central 
plane for a Rosseland optical depth of ~0.7 and the electron density (cm~^) at the same 
level. The column headed by log g is the log gravity (cgs units) in the z direction at a 
Rosseland optical depth of ~0.7. The column headed by tross is the Rosseland optical depth 
at the central plane. We call attention to the fact that the annuli are optically thick to the 
outer radius of the accretion disk. 
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Following; calculation of ann ulus models for assigned M^d and M, program SYNSPEC(v.48)3 
(iHubeny. Lanz. fc Jefferylll994l ) is used to produce a synthetic spectrum for each annulus. 
We adopted solar composition for all synthetic spectra. The synthetic spectra include con- 
tributions from the first 30 atomic species, assumed to be in LTE. 

We set SYNSPEC to produce light intensities at 10 zenith angles (equally spaced in 
cos(7) between 1.0 and 0.1, where 7 is the angle between the given direction and the local 
normal) at each calculated wavelength for all the individual annulus synthetic spectra. This 
choice automatically accounts for wavelength-dependent limb darkening in calculating syn- 
thetic spectra at different orbital inclinations. BINSYN calculates a synthetic spectrum for 
each of the annuli specified in the BINSYN system model by interpolation among the array 
of TLUSTY annulus models, with proper allowance for inclination. 

BINSYN models the complete CV system, including the WD, secondary star, accretion 
disk face, and accretion disk rim as separate entities. The model represents phase-dependent 
and inclination-dependent effects, including eclipse effects and irradiation effects, on all of 
the system objects. The representation of the stars requires polar Tefj values and gravity- 
darkening exponents. BINSYN represents the accretion disk by a specified number of annuli, 
where that number typically is larger than the number of TLUSTY annulus models. The 
accretion disk T^s profile may follow the standard model, but, alternatively, the profile may 
be specified by a separate input file. 

BINSYN has provision to repr e sent a bright spot on the rim face, but the low amplitude 



V light curve by iRibiero fc DiazI (120071 ) (Fig. 10) is consistent with light variation of an 
irradiated se condary star so we do not include a rim bright spot. Based on the analysis by 
SmakI (I2OO2I ) we set the rim Tes equal to the Tes of the outermost annulus. 



5. A V3885 Sgr system model 



We adopted system parameters listed in Table 4 for an initial test. Based on the simi- 
larity to IX Vel (comparable pe riods, visual magnitudes, parallaxes, CV class, spectroscopic 
properties (IHartley et al.l |2002| )). we first tried a standard rn odel accretion disk with the 
same M, 5.0x10"^ Mq yr'^ as for IX Vel Jlinnell et al.lboOTh . and, as for IX Vel, the (ini- 
tial) model included a 60,000K WD. In constructing a model accretion disk it is necessary 
to allow for the change in WD r a dius f rom the zero temperature model. We used Table 4a 
of IPanei. Althaus. fc Benvenutol J200oh to estimate the 60,000K WD radius at O.O134i?0. 



•^httpi/Znova. astro. umd.edu 
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Our BINSYN model specified 33 annuli for tlie accretion disk, witli assigned radii and cor- 
responding standard model Teg values listed in Table 6. The larger hot WD radius than 
a zero temperature WD radius at a given multiple of the 'WD radius' produces a lower 
Table 6 temperature than Table 5 (e.g., at the r/r„d value 6.0 in column 1 of both ta- 
bles). The synthetic spectrum for a given BINSYN annulus is calculated by interpolation 
(temperature-wise) among the Table 5 entries. Ideally, we could recalculate the whole series 
of annulus models (via TLUSTY) based on the new value of -Rwd, but the change in the 
individual annulus models would be very small and would not be warranted in view of other 
uncertainties (e.g., the system parallax). 

Table 5 ends at a Tefj of 7711K. TLUSTY fails for larger radii, even in the gray model 
case. The net effect is that the calculated BINSYN outer annuli have synthetic spectra 
corresponding to 7711K even if the specified standard model calls for a lower temperature. 
The effect on the synthetic spectrum of failing to follow the standard model is negligible. 
Adding a lower temperature (3500K) stellar synthetic spectrum to the array of Table 5, 
thereby permitting interpolation to a lower temperature in the outer BINSYN annuli, made 
no detectable change in the system synthetic spectrum. As se parate issues, the outer a nnuli 
temperatures should not fall below 6000K to avoid outbursts ( 10sakilll996l : lLasotall200ll ) and 
tidal dissipat ion is expect ed to raise the outer accretion disk temperature above the standard 
model value (ISmakll2002l ). 



The overall SED fit to the combined spectrum of §2.3, using initial parameters, and 
judged by visual estimate, was very good with a standard model accretion disk temperature 
profile. This was in contrast to IX Vel which required m odification of the te mperature pro- 
file from a standard model to achieve an acceptable fit (ILinnell et al.l 120071 ). Assessing the 
quahty of fit for V3885 Sgr is complicated by the presence of strong a bsorption features o f 
highly ionized material (Figure 1), probably associated with a wind (IHartley et al.ll2002l ). 
At this initial stage, the extreme synthetic spectrum FUV region had slightly too large flux: 
the only parameter (either M or Tefj^wd) available to adjust the flux while preserving the 
fixed normalizing factor was the WD Teg. (M could be adjusted but that route produced 
detectably larger residuals longward of 1500A. Small adjustments of M move the entire syn- 
thetic spectrum up or down without appreciable change in the spectral gradient. Similarly, 
a change in i does not change the spectral gradient. We consider a boundary layer in §6.) 
We reduced the WD T^e to 57,000K and found this produced an optimum fit as judged 
by eye estimate. The new model neglected the very small change in WD radius between 
Tgfj = 60, OOOK and 57,000K, which would have affected the tabular radii in Table 6. 



Figure 3 shows the FUV fit to the combined FU SE and STIS spectra, dereddened for 
E[B — V) = 0.02. We call attention to the broad absorption lines which we attribute to a 
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wind/chromosphere and that affect a major part of the FUSE spectrum. Note the deep 
ISM Lyman series hnes, the broad and deep absorption by 01 and NIII near 990A, the deep 
and possible mult i- component Ly/5 line near 1025A, and the broad lines of Silll, PV, SIV 
and SilV in the lllOA to 1140A interval. The overlapping and possibly multicomponent 
lines depress the apparent continuum over the entire interval from 970A to 1170A. 

The normalizing factor to divide the synthetic spectrum to superpose it on the observed, 
dereddened combined FUSE and STIS spectrum, corresponding to an exact distance of 110 
pc, is 1.1473x10"^^ cm^. We stress that this is not an adjustable parameter. In §5.1 we 
consider the effects of varying the 'exact' distance corresponding to ±la in the Hipparcos 
parallax. The spectrum at the bottom is the contribution of a 57,000K WD. For simplic- 
ity we assume the WD is rotating synchronously; we do not fit any observed absorption 
hnes in detail. The program assumes that the accretion disk blocks the view of the lower 
half of the WD. In addition, BINSYN allows for the partial visibility of that part of the 
WD not obscured by the accretion disk and above the observer's horizon, and includes 
wavelength-dependent limb darkening. The spectral resolution of all the synthetic spectra 
is O.O2A (§5.2.). The adopted orbital phase for the synthetic spectra was 0.75, in agree- 
ment with the observed spectra, §2. Adopting an orbital phase of 0.25 produced a clear 
misalignment of the calculated 960A peak with the observed one. 

We check the apparent best fit, by eye estimate, with a reduced chi squared (x^) 
test. We did not smooth the observed spectrum for this test. The presence of unmod- 
eled wind/chromosphere features in the observations produces large residuals (especially in 
the FUV, Figure 3). A test of our model fit to the observed continuum requires masking 
those spectral regions excluded from the fit (i.e., the wind/chromosphere features). We ini- 
tially masked spectral lines by hand, including the deep Lyman lines. (The Lyman lines in 
the model, mostly from the accretion disk, are washed out by Keplerian rotation.) This first 
stage masking produced a major reduction (from 154.39 to 7.38) in as compared with the 
unmasked spectrum. There remained a major "valley" in a residuals plot, described in the 
following paragraph, with the largest (negative) residuals near llOOA. Our previous discus- 
sion has identified this effect with overlapping absorption features from a wind/chromosphere 
which depresses the entire spectral interval over the approximate range 970A to 1170A. BIN- 
SYN has no present capability to model a wind/chromosphere and further progress requires 
masking the affected spectral interval. We therefore masked the entire region from 970A to 
II70A, producing a further large reduction in (from 7.38 to 2.65). 

Figure 4 shows plots of the residuals in the three cases. The top plot, displaced upward 
by 12.0 ordinate units, is the residuals for the unmasked case. The middle plot, displaced 
upward by 4.0 ordinate units, is the partially masked case. The lower (undisplaced) plot 
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is the final set of residuals adopted for this study. In each case a horizontal line marks 
the undisplaced reference. The change from the top to middle plot mostly results from 
masking of individual spectral features, such as the Lyman lines, the CIV A 1548 feature, 
etc. The STIS spectrum appears to have a defect beginning at 1710A, and the three STIS 
spectra show more variability longward of the Hell A 1640 line than in other spectral regions. 
We terminated our calculated fit just longward of the Hell A1640 line. The middle plot 
clearly shows the residual depression after individual spectral features have been masked 
and the bottom plot shows the overall improvement in residuals after masking the 970A to 
1170A region. The improvement is even more apparent on individual plots with extended 
ordinate axes. 

The lower plot still has fairly large outliers. Although further selective masking could be 
done, we consider it dangerous to go very far in suppressing observations that disagree with 
the model. Residual absorption features could be masked but there are positive outliers that 
then would bias the solution. We anticipate that it will be possible to locate an optimum 
model with the chosen masked set of observations. The final masked version of the observed 
combined spectra includes 4332 tabulated wavelengths and corresponding flux values. 

It is particularly unfortunate for this analysis that most of the FU SE spectrum must 

be masked, but it is unavoidable: the masked region clearly is strongly affected by the broad 
and deep absorption lines of high excitation species as described in §2.3. The high excitation 
involved excludes the lines from production by the accretion disk itself or the WD. Further, 
ignoring the high excitation aspect for the moment, it would not be possible to simulate the 
depressed 970A-1170A spectral interval by adopting a modified accretion disk T^s profile. 
Our experience with non-standard model accretion disks is that the calculated SED profile 
can have its slope changed over an extended spectral region but it is not possible to put a 
"hole" in the SED. 

We calculated values for a grid of 34 models in the M, Tegf^wd plane, centered on 
the best visual estimate model fit to the masked observations. Table 7 lists those values 
plus values for 4 additional i values (to be discussed separately). The aim is to provide 
a sufficiently extended grid that a subset of values define a contour that encloses 
the region of minimum x^, given the adopted values of other parameters. Calculation of 
a given model is a time-consuming process and a compromise is necessary between total 
computational time and density of grid points. Our computed grid is sparse but our judgment 
is that it meets our needs, given the uncertainties of separately adopted parameters. We 
have used the IDL routine CONTOUR to produce contour plots for the Table 7 i = 65° 
values. The M = 5.0x10"^ Mq yr"^ Teff,wd = 57, OOOK model has the smallest tabular 
value and we take that model to be our solution for V3885 Sgr. We discuss this choice and 



- 13 - 



the topic of grid density, including its adequacy, in more detail in §6. 

We have sampled the M, Teff^„d plane for the case M„d = O.SMq with two new 
models: M = 5.0xl0~9 Mq yi'^ and Teff,wd=57,000K, producing x^=3.06, and M = 
5.0x10^^ Mq yr~^ and Tcfr,wd=52,000K, producing x^=3.81. The corresponding val- 
ues from Table 7 are 2.65 and 3.39 respectively. A preferred value of M„d is not possible 
by analysis of the existing data. The SED accretion disk model is not a strong function 
of M^d for the change from 0.7Mq to O.SMq. Said differently. Figure 5 remains roughly 
constant, with a small systematic shift for changes in M„d of about 0.1 or less. 

An improv e d val ue of M„d may be possible with further optical region spectroscopy. 
Ribiero fc Diaj ( 2007 ) find that the Ha emission line arises from irradiaton of the secondary 



while the (doubled) A6678 Hel emission line associates with the accretion disk. If an accretion 
disk model that produces secondary component emission lines were available, a fit to the 
Hel line could constrain i. 

Table 7 includes 4 models for different assumed values of i and with other parameters 
the same as the preferred model. There is a very sharp minimum of at the adopted 
i = 65°. At this inclination there is an eclipse of the outer accretion disk at orbital phase 
0.0 with photometric effects much smaller than the predicted variation due to irradiation 
of the secondary. Increasing the inclination to 70° greatly deepens the depth of accretion 
disk eclipse and barely fails to eclipse the WD; it would produce light curve effects that 



are not observed (iRibiero fc DiazI 120071 . Figure 10). However the observed light variation is 
complex and may include a contribution from a rim bright spot. An inclination of 71° would 
eclipse the WD and produce a major light curve variation. It is not possible to use the 
i-dependent variation to determine i. The calculated system flux is a strong function of i 
because of the changing projection of the fiux-dominating accretion disk on the plane of the 
sky. The fit of the model spectrum to the observations requires a normalizing factor that 
is given by the Hipparcos parallax. The sharp minimum of for i = 65° results from the 
joint consistency of that parallax and the adopted i value. Actual determination of i would 
require observational effects intrinsic to the system, i.e., eclipses. 

Until a model that includes simulation of the wind/chromosphere is available it is doubt- 
ful that substantial further reduction in can be achieved. Based only on the x^ plot. Fig- 
ure 5, we estimate the parameter errors in M, Tes^wd, and i as 2.0x10"^° M© yr~^, 3000K, 
and, from Table 7, 2.0° respectively. However, shifts in M to match the uncertainty in the 
Hipparcos parallax (§5.1.) require a larger M error of 2.0x10^^ M© yr~^ which we adopt for 
the system solution. Further, in recognition of the sparse grid of models and their associated 
X^ values, we increase the estimated Teff,wd error to 5000K. 
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Our final system parameters are in Table 8. Where appropriate we have provided esti- 
mated error limits. Figure 6 shows the final model compared with the combined observations 
(not the masked observations). The blue line is the contribution of the 57,000K WD. The 
lower of the two black lines is the contribution of the accretion disk, and the upper black 
line is the system model. 



5.1. The effects of Hipparcos parallcLx uncertainty 

It is possible to test the effect of uncertainty in the Hipparcos parallax by changing 
the normalizing divisor for superposing the synthetic spectrum on the combined FUSE and 
STIS spectra. Changing the parallax by la to 7.16 mas changes the normalizing divisor 
from 1.147x10^^^ cm^ to 1.857x10'^^ cm^. The synthetic spectrum is displaced downward by 
approximately two fiux units on Figure 3 (the lOOOA peak moves from an ordinate value 
of 10. to 8.) Compensating for the change would require increasing M to perhaps 6.5 
or 7.0x10^^ Mq yr^^. Correspondingly, changing the parallax to 11.06 mas changes the 
normalizing divisor to 7.784x10^° cm^ and displaces the synthetic spectrum upward by about 
two fiux units, from a 1000 A peak of 10 to 12. Compensating would require decreasing M 
to perhaps 4.0 or 4.5xl0~^ Mq yr"^ 



5.2. Correction for the ISM 



Recent papers have i ncluded corrections for line absorption by the ISM (IGodon et al. 



20081 : iLinnell et al.ll2008al ). A custom spectral fitting package is used to estimate the tem- 
perature and density of the interstellar absorption lines of atomic and molecular hydrogen. 
The ISM model assumes that the temperature, bulk velocity, and turbulent velocity of the 
medium are the same for all atomic and molecular species, whereas the densities of atomic 
and molecular hydrogen, and the ratios of deuterium to hydrogen and metals (including he- 



lium) to hydrogen are adjust able parameters . The model uses atomic data of iMortonI (12000 



20031 ) and molecular data of lAbgralll (120001). The mo lecular hydrogen transmission values 
have been checked against those of iMcCandlisd (120031 ). The calculation requires fairly high 
spectral resolution and we have calculated all of the TLUSTY-SYNSPEC synthetic spectra 
at a resolution of 0.02A; we used the same resolution in calculating the system synthetic 
spectra. The ISM parameters are in Table 9. Figure 7 shows the corrected version of part 
of Figure 6. The ISM model produces a good representation of the observed deep cores of 
the Lyman series lines and the Lyman cutoff, as well as some lines in other spectral regions 
and clearly shows that the ISM produces the observed Lyman lines. The broad, high exci- 
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tation absorption regions from the accretion disk corona, chromosphere, or wind, identified 
in Table 2 and Table 3, remain as prominent residuals. As stated above, the Lyman lines 
were masked in the fit to the model and the spectral region shortward of the Lyman limit 
was masked so only the Lyman region continuum participated in the model solution. This 
was proper since a fit to the ISM does not affect the choice of the model parameters listed in 
Table 8. There are small positive residuals from the model near the short wavelength end. 
This could be interpreted as evidence for an unmodeled very high temperature contribution, 
discussed in the following section. 



6. Discussion 



In a well-known paper IWadd (119881 ) showed that neither steady state model accretion 
disks based on Planck functions or stellar model atmospheres could simultaneously fit the 
colors and absolute luminosities of a set of NL systems ('V3885 Sgr was included in the 
study). Tomographic analyses of accretion disks (IRutten. van Paradijs. fc Tinbergenlll992l ) 
show clear departures from standard model temperature profile s (FKR). Recent analyses 
of the NL systems MV Lyr, IX Vel, QU Car, and UX UMa jLinnell et all hood . bopj . 
200830) showed, in each case, that a standard model accretion disk synthetic spectrum 
could not accurately fit observed spectra. These results suggest that physical conditions in 
possibly all NL systems are so complex, with winds, coronae, iron curtains, etc., that an 
accurate fit by a standard model synthetic spectrum is too much to expect. Consequently, 
the most interesting result of this study is that a standard model synthetic spectrum plus WD 
accurately fits the observations (x^=2.65) after masking the effects of a wind/chromosphere. 
This conclusion is particularly true in light of the fact that the system distance is known 
so the fit is in absolute fiux units. It is worth emphasizing that the model presented here 
includes a measurable contribution from a hot source, in addition to the accretion disk, 
sufficiently large that limits can be set on the hot source T^s. Assuming that the hot source 
is a WD (see below), this is the first time the effective temperature of a hot WD has been 
constrained during the high brightness state of a nova-like variable. 

As discussed by FKR, half of the potential energy liberated in the fall of the mass transfer 
stream from the LI point (essentially from infinity) to the WD surface appears as radiated en- 
ergy. Energy conservation requires an accounting for the other half and the presence of a hot 
boundary layer (BL) is a co mmon prescrip tion. The problem is that BLs with the prescribed 
properties are not observed (IC6rdovalll995l ). The abse nce of an observed but predicted (FKR) 



BL in high M CV systems h a s an extensive history : ( Ferland 



Patterson fc Ravmondlll985l : iHoare fc Drewlll99ll . Il993l : IVrtilek et al.lll994l : lldan fc Shaviv 



1982 



Kallman &: Jensen 



1985 
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19961 ). One proposed explanation is that the WD is rotating close to Keplerian rot ation, 
with no BL predicted. However, HST measurements starting with ISion et al.l (Il994j ) have 
measured rotational WD velocities too low to explain "missing" BLs. An alternative ex- 
planatory theme is that the a bsence of an observed BL associates with a wind, and in V3885 
Sgr a wind is clearly present (IHartley et al.ll2002l . |2005| ). 



BINSYN has an option to simulate a BL by specifying the innermos t annulus radial 
thickness and T^s. We applied the theory of lGodon. Regev. fc Shavivl (119951 ) to obtain those 
two parameters. This theory derives a BL with a Tegf~100, OOOK and rad ial thickness of 
abou t 0.1i?wd- This BL T^ff is consisten t with a val ue from wind ionization (IHoare &: Drew 
19911 ) and the theory of lldan fc Shavivl (119961 ). The iGodon. Regev. fc Shavivl (119951 ) model 
adopted M„d = I.OMq and Twd,eff=20,000K. Addition of the BL made a nearly negligible 
change in our system SED. There was a barely detectable increase in the system flux at the 
short wavelength limit over the SED of our original model, producing a slightly better fit 
in that spectral region. The evidence in favor of a BL is inconclusive; a model including a 
100, OO OK BL and a 57, OOOK WD is consistent with the observations. iPopham &: Narayan 
(I1995I ) calculate BL models and illustrate their theory for the case M = 10~^ Mq yr~^ 
and zero WD rotation, leading to a somewhat higher BL Teg. In a discussion among the 
coauthors Godon noted that the BL Teg might be as high as 200, OOOK but probably would not 
be visible at i = 65° because of an inflated inner disk; an additional reason for not including 
a BL. A model including the 100, OOOK BL but with our WD contribution suppressed (i.e., 
requiring the hot source to be the BL exclusively) would be strongly inconsistent with the 
observations. We are unaware of any BL models that treat MHD effects. 

Is it possible to exclude a different source for the WD contribution of our present model? 
We seek some spectral feature that could distinguish the WD contribution from a different 
prospective source. In Figure 6 the WD Ly 7 absorption line appears to make an important 
contribution to the system synthetic spectrum. A different smoothly varying source, such 
as a black body (BB) or power law source might produce a much poorer fit to the observed 
Ly 7 wings and so be excludable. We tested this possibility by replacing the WD synthetic 
spectrum by a 57,000K BB with the same size as the replaced WD. The existing model 
uses light intensities at 10 zenith angles (§4), produced by SYNSPEC, for each annulus. 
Use of the BB spectrum required a change to fiux-based synthetic spectra with an adopted 
wavelength- independent limb darkening coefficient; we chose 0.6. The system synthetic 
spectrum SED was essentially unchanged but required a slight adjustment of the normalizing 
factor to reproduce the same fit as with the original model, likely a result of the change to 
a wavelength- independent limb darkening coefficient. The fit at the short wavelength limit 
was slightly better, but it is uncertain whether this arises from the BB replacement or the 
change in limb darkening. In any event, the Ly 7 line was slightly shallower in the new 
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simulation, arguing weakly in favor of the presence of a WD, but the fit to the wings was 
essentially the same. The proposed spectral feature test fails and we are unable to exclude a 
BB as the source of the required hot addition to the accretion disk. The SED of a power law 
source would differ from that of a BB source and should produce a poorer fit to the observed 
spectrum; without explicit test we discount this type source as a reasonable candidate. 

To answer the question "What is the source of the hot addition to the accretion disk?" 
the most reasonable working hypothesis is that the hot source is a WD. A WD is known to 
be present (it is a CV system) and the inclination, i, is such that the WD should be visible. 
A spectral source substitute for a WD (e.g., BB) would need not only to explain why the 
WD spectrum is not seen but also why a BB source of the same Teg and size (which sets the 
radii of the accretion disk annuli and thus their Tgfj values (Table 6)) and physical location 
(exactly centered on the accretion disk) is a reasonable substitute. 

The adequacy of our set of 34 models in populating the x^(M, Tgfr^wd) plane requires 
justification At issue is the question "How dense a grid is adequate?" Phrased differently, 
"As we produce a denser and denser grid, and keeping in mind the committment of computer 
time, how do we know when to stop?" A very important point is that our solution adopts 
values for the parameters M(wd), M(sec), Period, i, parallax, and E{B — V). A solution with 
a different value of any one of these would require recalculation of an entire new plane of 
X^(M, Teff^wd) values, like Figure 5, and search of it for the minimum x^. That process would 
be required, for example, if we study the implications of changing the adopted parallax by 
one standard error. The solution we have found from our set of values is valid only under 
the assumption that the other fixed parameter values are accurate. 

Our Table 7 values for a specific M and variable Teg^wd, or specific Teg^wd and variable 
M, extend to a x^, when rounded to the nearest half integer, of at least 4.0. We consider 
that value to be enough larger than the minimum of 2.65, and the resolution in M and Tgfj ^d 
sufficiently fine, to provide reasonable assurance that the (ideal) global minimum has been 
encircled. Figure 5, produced with IDL routine CONTOUR, presents contours of constant 
including closed contours that encircle the tabular minimum. 

We reiterate the point that possibly our most interesting result is the accurate fit to a 
standard model. We note that a visual estimate (" chi by eye" ) selected the best model. 



Hartley et al.l (120021 ) discuss the wind in V3885 Sgr and point to the lack of correla- 
tion between wind activity and system luminosity, an inconsistency with expectation for 
line-driven winds. The theory of line- driven winds h as been applied to CV systems by 



Perevra. Hillier. &: Turnshekl (120061 ) and iProga J2003al lb[). See the review by iDrew &: Progal 



(I2OOOI ). iProgal (l2003al ) considers both radiation-driven and MHD effects and shows that hy- 
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brid models that i nclude both M HD and line-driving produce larger Mwind values, of order 
IQ-^^ Mq yi-\ but IProgal fl2003bf ) shows that a pure lin e- driven wind p roduc es an improved 
fit to observed line shapes but with a lower Mwind- As [Hartley et al.l (120021 ) point out, the 
lack of correlation between wind activity and system luminosity may argue in favor of MHD 
effects. In either case, the theoretically derived values of M^md are far smaller than the mass 
transfer M and essentially all of the mass transfer stream must be deposited on the WD. 



Using eq. 2 of ISionI (Il995l ). the Teg of a compressionally- heated WD is given by T^fi = 
[(/GM„dM)/(47ri?^j(j)]°'^^, where / is the fraction of the accretion energy applied to com- 
pressional heating and is in the range 0.15-0.20. For / = 0.15 and the parameters of Table 8 
we calculate Tefr=52,500K, while for / = 0.20 we find Tefr=56,414K, well within the estimated 
error limit of our 57,000K WD. 

This study raises an obvious question: Why does a standard model accretion disk repre- 
sent V3885 Sgr accurately while other apparently comparable systems violate the standard 
model? To rephrase the question: What equation or equations defining the standard model 
do the nonconforming systems violate? The standard model equations (e.g., FKR eq. 5.68) 
adopt axial symrnetry a nd hydrostatic equilibrium, avoid MHD considerations (the MRI 
(iBalbus fc Hawleyl Il99ll ) source of viscosity), and assume a constant a for the entire ac- 
cretion disk, corresponding to a const ant local M within the accretion disk. The observed 
presence of spiral waves in V3885 Sgr (IHartley et al.ll2005l ) violates the axial symmetry con- 
dition even in a system that fits the standard model. The impact of the mass transfer stream 
produces a departure from axial symmetry. These phenomena illustrate the point that the 
standard model is a simplified but very useful idealization. Just as weather in the Earth's 
atmosphere, from physically complex phenomena, produces transient but extended local de- 
partures from the ideally accurate hydrostatic equilibrium, we speculate that transients of 
unknown but physically complex origin may occur in localized regions of accretion disks 
that perturb the radiation characteristics from the standard model average. lUE spectra 
of UX Ursae Mai ori s show substantial time- dependent SED changes for that NL system 
(ILinnell et al.l l2008bl ). The adherence of V3885 Sgr to the standard model in this study 
may indicate an infrequently violated average. On the other hand, the observation times 
of nonconforming systems may occur at times when transients are prevalent but detailed 
frequent monitoring might discover times when their SEDs fit the standard model. 



7. Summary 



We use the iKniggd (120071 ) period-secondary mass relation to assign the secondary mass 
and the mean WD mass from the same paper, consistent with an observationally-determined 
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mass ratio, to set the V3885 Sgr WD mass. Based on a analysis using 38 system models 
and corresponding synthetic spectra, we find that a standard model accretion disk with 
M = 5.0±2.0xl0-9 Mq yr-^ a 57, 000±5000K 0.7Mq WD, and an orbital inclination of 
65±2° provide a system synthetic spectrum which optimally fits observed FU SE and STIS 
spectra of V3885 Sgr. The result is robust for small changes in the assumed M^d- A 
Hipparcos parallax fixes the divisor enabling the system synthetic spectrum to be matched 
to the observed spectra and establishes the fit on an absolute flux basis. The model thus 
provides realistic constraints on M and T^s for a large M system above the period gap. 
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Table 1: FUV Observations of V3885 Sgr: HST fc FUSE Spectra 



Date 


Telescope 


Exp time 


Dataset 


Filter/Grating 


Operation 


Calibration 


(dd/mm/yyyy) 


/Instrument 


(sec) 




/ Aperture 


Mode 


Software 


24 May 2000 


FUSE 


12392 


P1870101000 


LWRS 


TTAG 


CalFUSE 3.0. 


30 Apr 2000 


STIS 


2195 


O5BI04010 


E140M/0.2x0.2 


TIME-TAG 


CALSTIS 2.2 
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Table 2. FUSE lines 



Line 


Wavelength 


Ao 


Origin 


Identification 


(A) 


(A) 




H T 

11 1 


Q14 1 


Q14 9Q 


IGTTl 


Hi 


Q14 43 


914 58 


1*5 TTl 


Ht 


Q14 7Q 


914 92 


lolll 


Ht 


Q15 1Q 


91 33 

ly XtJ . 


lolll 


H T 

11 i 


Q1 fiS 


QI ^ 89 


loili 


Ht 




916 43 

t/ X\J 


lolll 


O T 




916 82 


lolll 


H T 

11 1 


Q17 n"! 


917 18 


lolll 


H T 

11 1 


Q18 nn 

zf ±o .yjyj 


918 1 3 

c/±0. xO 


lolll 


Ht 


Q1Q 22 


919 35 


ICTTl 
lOlll 


Ht 


Q2n 84 


920 96 


lolll 


H T 

11 i 




991 86 

c/Zj X .OU 


lolll 


Ht 


Q2S ni 


993 1 5 


lolll 


N TV 
ly 1 V 




993 9 




n T 


Q94 8S 


994 95 


lolll 


O T 




925 45 


lolll 


H T 

11 1 


Q2fi 1 1 


926 23 


lolll 


Ot 


Q2Q 41 


929 52 


lOlll 


Ot 




930 26 


ICTTl 
lOlll 


lU T 

n I 


ydu.oi 


non Tc: 
9oU. / O 


ism 


Svi 


~932 


933.38 




Hi 


936.51 


936.75 


ism 


Hi 


937.67 


937.80 


ism 


Svi 


~943 


944.52 




Oi 


948.57 


948.69 


ism 


Hi 


949.61 


949.74 


ism 


Oi 


950.76 


950.88 


ism 


Ni+Oi 


952.26 




ism 


Nl 


953.30 


953.42 


ism 
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Table 2 — Continued 



Line Wavelength Aq Origin 

Identification (A) (A) 



Nl 


953.54 


953.65 


ism 


Nl 


953.87 


953.97 


ism 


Nl 


953.97 


954.10 


ism 


Pll 


962.44 


962.57 


ism 


Pll 


963.68 


963.80 


ism 


Nl 


963.87 


963.99 


ism 


Nl 


964.50 


964.63 


ism 


Nl 


964.94 


965.04 


ism 


Oi 


971.61 


971.74 


ism 


Hi 


972.34 


972.54 


ism 


Cm 


~976 


977.02 


s 


Cm 


976.88 


977.02 


ism 


Oi 


976.32 


976.45 


ism 


Oi 


988.56 


988.65 


ism 


Sin 


989.75 


989.87 


ism 


Nm 


~990 


991.5 


s 


Sim 


997.26 


997.38 


ism 


Sin 


1020.63 


1020.70 


ism 


Hi 


1025.37 


1025.72 


s,ism 


Hi 


1025.68 


1025.72 


s,ism 


Ovi 


1030.7 


1031.93 


s 


Ovi 


1036.4 


1037.62 


s 


Cii 


1036.24 


1036.34 


ism 


Oi 


1039.18 


1039.23 


ism 


An 


1048.15 


1048.22 


ism 


Siv 


1062.3 


1062.65 


s 


An 


1066.59 


1066.66 


ism 


Siv 


1072.5 


1072.97 


s 


Nil 


1083.74 


1083.99 


ism 
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Table 2 — Continued 



Line 


Wavelength 


Ao 


Origin 


Identification 


(A) 


(A) 




Sim 


1109.4 


1108.36 


s 






1109.94 


s 


Sim 


1113.0 


1113.2 


s 


Siv 


1117.3 


1117.76 


s 


Pv 


1117.3 


1117.98 


s 




1127.4 


1128.01 


s 


Si IV 


1122.0 


1122.5 


s 




1127.4 


1128.3 


s 


Ni 


1134.08 


1134.16 


ism 


Nl 


1134.33 


1134.42 


ism 


Nl 


1134.89 


1134.98 


ism 


Cm 


1175.0 


1174.9-1176.4 


s 



Note. — a=possibly with terrestrial contamination; 
ism=ISM; s=source 
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Table 3. STIS lines 



Line Wavelength Aq Origin 

Identification (A) (A) 



Cm 


1173.62 


1174.9-1176.4 


s 


Sin 


1190.36 


1190.42 


ism 




1193.24 


1193.29 


ism 


Ni 


1199.49 


1199.55 


ism 




1200.18 


1200.22 


ism 




1200.67 


1200.71 


ism 


Sim 


1204.77 


1206.50 


s 


Sim 


1206.45 


1206.50 


ism 


Hi 


1215.6 


1215.67 


s 


Nv 


1236.9 


1238.82 


s 




1240.3 


1242.80 


s 


Nv 


1238.63 


1238.82 


ism 




1242.59 


1242.80 


ism 


Sii 


1250.55 


1250.58 


ism 




1253.77 


1253.81 


ism 




1259.45 


1259.52 


ism 


Sin 


1260.37 


1260.42 


ism 


Sim 


1297.3 


1294.55 


s 






1296.73 


s 






1298.89 


s 


Oi 


1302.12 


1302.17 


ism 


Sin 


1304.33 


1304.37 


ism 


Cii 


1334.47 


1334.53 


ism 




1334.47 


1334.66 


ism 




1335.66 


1335.71 


ism 


Si IV 


1391.4 


1393.76 


s 




1400.35 


1402.77 


s 


Sin 


1526.62 


1526.71 


ism 


Civ 


~1545 


1548.20 


s 
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Table 3 — Continued 



Line 


Wavelength 


Ao 


Origin 


Identification 


(A) 


(A) 








1550.78 


s 


Civ 


1547.95 


1548.20 


ism 




1550.54 


1550.78 


ism 


Hen 


1637.7 


1640.5 


s 


Alii 


1670.73 


1670.8 


ism 



Note. — ism=ISM; s=source; there is possi- 
bly some atmospheric absorption contamination 
for Ni, Nil, Hi, Oi. 
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Table 4. V3885 Sgr Initial System Parameters 



parameter value parameter value 



M^d O.7OM0* i 65° 

M2 0A75Mq^ M S.OxlO-^Meyr-i " 

P 0.20716071 day ^ d 110±30 pc ^ 



Kniggel (120071) mean M..^: iHartlev et all (120051 ) . 0.55Mq to 



0.8M^ 



0) 



Ribiero fc DiazI (120071 ) . 0.3Mq to 0.9M£ 







Kniggd (I2OO7I ) cali brated P : M 2 rela tion; [Hartley et al. 



(120051 ). g=0.6 to 0.8; iRibiero fc DiazI (120071 ). M2 = O.2M0 to 
O.55M0 



Ribiero fc DiazI (120071 ) 



iHaug fc Dreschell (119851 ). i between 60° and 70 
Ribiero fc DiazI J2007[ l. i between 45° and 75°; Iflartlev et at 



(I2OO5I ). i > 65"= 



Linnell et al.l (120071 ) and similarity to IX Vel (iHartley et al. 

fcooj) 

^Hipparcos value 
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Table 5. Properties of accretion disk with mass transfer rate M = 5.0x10 ^ MqJt ^ and 

WD mass of 0.70Mc.. 









mo 


log 


9 










^Ross 


1.36 


53136 


1 


221E4 


6 


84 


8 


27E7 


2 


0E17 


Z 


ooJi/4 


2.00 


47642 


1 


491E4 


6 


58 


1 


43E8 


1 


3E17 


2 


78E4 


3.00 


38525 


1 


456E4 


6 


28 


2 


43E8 


7 


7E16 


3 


35E4 


4.00 


32380 


1 


356E4 


6 


06 


3 


49E8 


5 


1E16 


3 


92E4 


5.00 


28086 


1 


259E4 


5 


89 


4 


60E8 


3 


8E16 


4 


52E4 


6.00 


24918 


1 


176E4 


5 


75 


5 


76E8 


3 


0E16 


5 


16E4 


8.00 


20531 


1 


043E4 


5 


53 


8 


14E8 


2 


5E16 


6 


56E4 


10.00 


17612 


9 


439E3 


5 


35 


1 


06E9 


1 


9E16 


8 


05E4 


12.00 


15512 


8 


664E3 


5 


21 


1 


32E9 


1 


5E16 


9 


53E4 


14.00 


13923 


8 


040E3 


5 


09 


1 


58E9 


1 


2E16 


1 


08E5 


16.00 


12670 


7 


526E3 


4 


98 


1 


85E9 


9 


9E15 


1 


18E5 


18.00 


11653 


7 


092E3 


4 


89 


2 


12E9 


8 


3E15 


1 


23E5 


20.00 


10810 


6 


721E3 


4 


80 


2 


38E9 


6 


7E15 


1 


23E5 


24.00 


9487 


6 


115E3 


4 


65 


2 


92E9 


5 


1E15 


1 


17E5 


28.00 


8491 


5 


638E3 


4 


51 


3 


57E9 


4 


7E15 


1 


07E5 


32.00 


7711 


5 


251E3 


4 


46 


4 


90E9 


2 


2E14 


9 


32E4 



Note. — Each line in the t able represents a separate annulus. 
A IShakura fc SunyaevI (119731 ) viscosity parameter a = 0.1 was 
used in calculating all annuli. The WD radius, r^d,o, is the ra- 
dius, O.OlO91i?0, of a zero temperature Hamada-Salpeter carbon 
model. See the text (54.) for a discussion of the table units. 
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Table 6. Temperature profile for V3885 Sgr accretion disk with mass transfer rate of 
M = 5.0x10-9 MqJy-^ and WD mass of O.TOM©. 











r/r^d 


Teff 


1.00 


43608 


13.32 


12308 


27 


94 


7260 


1.18 


43608 


14.65 


11511 


29 


27 


7021 


1.36 


45358 


15.98 


10826 


30 


59 


6799 


2.69 


34985 


17.31 


10229 


31 


92 


6594 


4.02 


27556 


18.64 


9705 


33 


25 


6402 


5.35 


22948 


19.96 


9240 


34 


58 


6223 


6.68 


19807 


21.29 


8824 


35 


91 


6055 


8.01 


17518 


22.62 


8449 


37 


24 


5898 


9.33 


15767 


23.95 


8110 


38 


57 


5749 


10.66 


14380 


25.28 


7802 


39 


90 


5610 


11.99 


13249 


26.61 


7520 


41 


22 


5478 



Note. — The WD radius r^d is the radius of 
a 60, 000K 0.70M^ WD interpolate d from Table 
4a of iPanei. Althaus. fc Benvenutd (120001 ). The 
Teff = 43608K for r/rwd=l-00 refers to the inner 
edge of the innermost annulus. The same Teff for 
r/r„d= 1-18 refers to the outer edge of the inner- 
most annulus which coincides with the inner edge 
of the next annulus. The remaining Teff values re- 
fer to the inner edge of the corresponding annuli. 
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Table 7. Values of as function of model parameters 



M (xlO ^Mq/jt:) 




ndeg.) 


~2 

X 


M (xlO ^Mq/yt) 


rwd(K) 


ndeg.) 


~2 


A n 


7(1 nnn 


65 


7 51 


5 

\j .yj 


57 nnn 


64 


3 Q5 


4 


75 (1(1(1 


65 


5 60 

•\J\j 


5 


57 nnn 


66 


4 43 




57 (1(1(1 


65 


6 86 


5 


57 nnn 


70 


3Q 73 


4.5 


60,000 


65 


5.05 


5.2 


40,000 


65 


4.69 


4.5 


65,000 


65 


3.93 


5.2 


45,000 


65 


3.14 


4.5 


70,000 


65 


3.14 


5.2 


50,000 


65 


2.91 


4.5 


75,000 


65 


3.97 


5.2 


52,000 


65 


2.73 


4.8 


52,000 


65 


5.22 


5.2 


57,000 


65 


3.n4 


4.8 


57,000 


65 


3.42 


5.2 


6n,nnn 


65 


4.41 


4.8 


60,000 


65 


3.00 


5.5 


35,nnn 


65 


4.34 


4.8 


65,000 


65 


3.40 


5.5 


40,000 


65 


3.08 


4.8 


70,000 


65 


4.26 


5.5 


45,000 


65 


3.19 


5.0 


45,000 


65 


4.54 


5.5 


50,000 


65 


3.24 


5.0 


50,000 


65 


4.25 


5.5 


52,000 


65 


3.93 


5.0 


52,000 


65 


3.39 


5.5 


57,000 


65 


5.93 


5.0 


57,000 


65 


2.65 


6.0 


35,000 


65 


4.03 


5.0 


60,000 


65 


3.44 


6.0 


40,000 


65 


6.48 


5.0 


65,000 


65 


5.80 


6.0 


45,nnn 


65 


4.13 


5.0 


57,000 


60 


49.94 


6.0 


50,000 


65 


9.88 
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Table 8. V3885 Sgr Model System Parameters 



parameter 


value 


parameter 


value 


71 /T 


(J.7U±U.1M0 


rri / 1 \ 

Tcfr,s(pole) 


oo50±100K 


71 T 

Ms 


O.475M0 


Teff,s(pomt) 


2214K 


M 


5.O±2.OxlO-9M0yr-i 


Teff,s(side) 


3598K 


P 


0.2071607 days 


Teff,s(back) 


3489K 


D 


1.55396i?0 


rs(pole) 


O.2O9i?0 


^wd 


116.0 


(point) 


O.296i?0 




3.2053112 


rs(side) 


O.218i?0 


i 


65±2° 


r^back) 


O.239i?0 


^cflf.wd 


57,000±5000K 


log fi's(pole) 


4.73 


'^wd 


O.O134i?0 


log (point) 


-.390 


log g„d 


8.02 


log (side) 


4.65 


^wd 


1.0 


log g's(back) 


4.48 


As 


0.6 




O.56i?0 


/3wd 


0.25 




O.OlO91i?0 




0.08 


H 


O.lOi?0 



Note. — wd refers to the WD; s refers to the secondary star. 
D is the component separation of centers, 1] is a Roche potential. 
Temperatures are polar values, A values are bolometric albedos, 
and P values are gravity- darkening exponents. Va specifies the outer 
radius of the accretion disk, set at the tidal cut-off radius, and rb 
is the accretion disk inner radius, as determined in the final system 
model, while H is the semi-height of the accretion disk rim. 



Table 9. Parameters of ISM model 



Parameter Value 



N(H) l.OxlOi^cm-2 

D/H 2x10-^ 

N(H2) 1x10^^-2 

vel. -50.0 km s"^ 

temp. 250 K 

metallicity 1.0 (1.0 — solar) 
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Wavelength (A) 

Fig. 1. — FUSE plus ST IS combined spectra with line identifications. 
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1140 
12 



B 

o 

'(/} 
lao 

CD 



X 
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1145 1150 1155 1160 

Wavelength (A) 



1165 



1170 




1000 



1600 



1200 1400 
Wavelength (A) 

Fig. 2. — Combined FUSE and STIS spectra. The upper panel shows the fit of the spectra 
in the overlap region when small normalizing factors are applied to the STIS spectra. The 
FUSE spectrum is red; the first STIS spectrum (Table 1), divided by 0.98, is cyan; the 
second STIS spectrum, divided by 1.12, is magenta; the third STIS spectrum, divided by 
0.94, is brown. Note the accurate superposition of the STIS spectra with application of 
normalizing factors. 
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900 1000 1100 1200 1300 

Wavelength (A) 

Fig. 3. — FUV fit of system synthetic spectrum to FUSE and STIS combined spectrum 
(gray plot) of V3885Sgr. Note the extensive absorption by local highly ionized material. 
The lower spectrum is the contribution of a 57,000K WD; the top synthetic spectrum is the 
system synthetic spectrum including the WD and a M — 5.0x10"^ Mq yr~^ standard model 
accretion disk. See the text for a discussion. 




1000 1200 1400 1600 

Wavelength (A) 

Fig. 4. — Plots of residuals from the visually selected best model. The top plot, displaced 
upward by 12.0 ordinate units, is the residuals in the unmasked case. The middle plot, 
displaced upward by 4.0 ordinate units, is the residuals for a partially masked observed 
spectrum. The bottom plot is the residuals for the finally adopted masked case. Note the 
three horizontal lines marking residual values of 0.0 for the three cases. See the text for a 
discussion. 
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Fig. 5. — Contour plot showing values of in the M, Tes,wd plane for the system models 
tabulated in Table 7. An asterisk marks the value for the preferred solution. See the text 
for a discussion. 
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Fig. 6. — As in Figure 3 but for the full wavelength range of the combined spectra. The 
synthetic spectrum of the accretion disk alone is the next-to-top synthetic spectrum. The 
uppermost system synthetic spectrum is the sum of the accretion disk and the WD (blue 
line). The contributions of the secondary star and accretion disk rim are negligible. 
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Fig. 7— The model of Figure 6 but with ISM correction. Note the cutoff at 912A. 



